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, The asymmetrically mounted flat plasma actuator is studied using a self-consistent 2-DIM fluid 

£NJ ■ model at atmospheric pressure. The computational model use the drift-diffusion approximation and 

a simple plasma phenomenological kinetic model. It is investigated its electrical and kinetic proper- 
ties, and calculated the charged species concentrations, surface charge density, electrohydrodynamic 
forces and gas speed. The present computational model contributes to understand the main physical 
mechanisms and methods to improve its performance. 
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C3 i I. INTRODUCTION 

There has been a growing interest in the field of plasma aerodynamics related to its outstanding importance in active 
flow control, overriding the use of mechanical flaps [H, [2, H, 0| • Plasma actuators create a plasma above a blunt body 



that modify the laminar-turbulent transition inside the boundary layer [H, Q , even at a high angle of attack 0, 0] they 
induce or reduce the fluid separation, and thus reducing drag [l[ and increasing lift [1, G3- They also allow sonic 
boom minimization schemes [Til 12\, avoiding unwanted vibrations or noise (lS. fl4l|. sterilizing or decontaminating 
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Oh. surfaces and even frost removal [la ] on any airfoil, jet engine or wind turbine, through pure electromagnetic control. 
Its promising potential extend to flow control at hypersonic speeds flU [l7j (while still using a jet-reaction aircraft 
i , propeller). 

■ The asymmetric dielectric barrier discharge (DBD) plasma actuator is a normal glow discharge that, like all normal 
\ glow discharges, operates at the Stoletow point. This guarantees that the generation of the ion-electron pairs at one 
atmosphere is done efficiently. In air, the minimum energy cost is 81 electron volts per ion-electron pair formed in 
the plasma. In plasma torches, this energy cost can be of the order of one keV/ion-electron pair; in arcs, it can range 
. . from 10 to 50 keV/ion-electron pair. 

Particle and fluid simulations have been done [HI, [l!| [2(| for a plasma actuator in pure oxygen and pure nitrogen 
showing the formation of an asymmetrical force that accelerates the ions dragging the neutral fluid in the direction 
of the buried electrode. A net force arises because the plasma density, and consequently, momentum transfer are 
greater during the second half of the bias cycle, due partially to the ion density greater a factor of 10 times during 
the second half-cycle. Thus, in each cycle there is created a total unidirectional force towards the buried electrode 
that can create neutral fluid flow velocities on the order of 8 m/s. 

Two dimensional fluid models of a DBD plasma actuator have been made [U, [52] which calculate the total force 
on ions and neutral particles and show that the force generated is of the same nature as the electric wind in a corona 
discharge, with the difference that the force in the DBD is localized in the cathode sheath region of the discharge and 
expands along the dielectric surface. While the intensity of this force is much larger than the existing force of a dc 
corona discharge, it is active during less than a hundred nanoseconds for each discharge pulse and, consequently, the 
time averaged forces are of the same magnitude in both cases. 

The use of voltage pulses in plasma actuators, by modulation of the high frequency excitation voltage carrier wave 
by a square wave, introduces mean and unsteady velocity components and thus the air momentum is composed of 
both time- mean and oscillatory components of momentum, the overall effect improving momentum efficiency [23j | . 
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FIG. 1: Schematic of the asymmetric plasma actuator. 

Actuators placed on the leading edge of an airfoil can control the boundary layer separation, while if located at the 
trailing edge can control lift [2(| ■ Enloe et al. |24l [25| have found experimentally that the thrust T and maximum 
induced speed u max are proportional to the input power P, which depends nonlinearly on the voltage drop AV across 
the dielectric T oc u max ocPa AV 7 ' 2 . 

The aim of this paper is to present a self-consistent two dimensional modeling of the temporal and spatial devel- 
opment of asymmetric DBD plasma actuator using an EHD code (CODEHD) developed in our group. The com- 
putational model solves the governing equation in the drift-diffusion approximation and a plasma phenomenological 
kinetic model. Fig. [1] shows the rather simple configuration of the plasma actuator in coplanar configuration. 

II. NUMERICAL MODEL 
A. Description 

At our knowledge the first comprehensive kinetic model of a dielectric barrier discharge plasma actuator was 
published by Singh et al. [26|. Gadri has shown that an atmospheric glow discharge is characterized by 
the same phenomenology as low-pressure dc glow discharge. To subdue numerical complexity no detailed plasma 
chemistry with neutral heavy species is presently addressed. At this stage it is only considered the kinetics involving 
electrically charged species supposedly playing a determinant role at atmospheric pressure: Nj, , 0^", O^, and 
electrons. From the charged species populations and as well the electric field controlling their dynamics, it can be 
studied electrohydrodynamics (EHD) with interest to plasma actuators, like the body forces acting on the plasma 
horizontally (neutral flow control) and perpendicularly (boundary-layer control) to the energized electrode; and the 
neutral particles average speed. 

The applied voltage has a sinusoidal wave form V(t) = Vd c + V rms sm(wt)/\/2, where the root mean square voltage, 
V rmsi in this case of study is 5 kV and the applied frequency is / = 5 kHz. Therefore, the dynamical time is T = 200 
/is. 

The plasma actuator simulation domain is a 2-Dim Cartesian geometry with the total length along the Ox-axis 
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TABLE I: List of reactions taken into account with CODEHD. Rate coefficients for chemical processes were taken from Ref. [2S 

Kind of reaction Process Rate coefficient 

Ionization e + N 2 -> 2e + N+ v™ 2 a 



Ionization e + 2 —> 2e + 2 v. 



lot: 
2 
ion 



3-body electron attachment e + 2 + 2 -> + 2 #ai = 1-4 x 1CT 29 (^) exp(-600/T 9 ).Ki(T fl , T e ) (cm 6 /s) 6 
3-body electron attachment e + 2 + N 2 ^ 2 + N 2 K a2 = 1.07 x 1()- 31 (^) 2 K 2 {T g , T e ) (cm 6 /s) c 

1-10. r 



Collisional detachment 2 + 2 -> e + 20 2 ^da = 2.9 x lO" 10 ^ ^ exp(-5590/T g ) (cm d /s) 

e-ion dissociative recombination TV^ + e — > 2A^ j3=2.8x 10~ 7 y^|p (cm 3 /s) 

e-ion dissociative recombination 2 + e — > 20 /3 = 2.8 x 10 _7 ^/lp2 (cm 3 /s) 

2-body ion-ion recombination + iV 2 + -> 2 + N 2 j3 u = 2 X 10" 7 ,/|^[1 + lO" 19 ^^) 2 ] (cm 3 /s) 

2-body ion-ion recombination O^T + C# -> 20 2 ft, = 2 x 10" 7 ./|^[1 + 10" 19 iV(^) 2 ] (cm 3 /s) 

Ion-conversion iV 2 + + iV 2 + 7V 2 -> iV^" + e # icl = 5 x 10 -29 (cm 6 /s) 

Recombination " N+ + e ^ 2N 2 K r2 = 2.3 x 10- 6 /(Te/300) - 56 

Ion-conversion N+ + N 2 -> 7V+ + 2iV 2 X lc2 = 2.1 x 10~ 16 exp(T s /121) (cm 3 /s) 

"Numerical data obtained by solving the quasi-stationary, homogeneous electron Boltzmann equation. See Ref. [29l | for details. 

6 With K x = exp(700(T e - T g )/(T e T g )). 

c With K 2 = exp(-70/T s ) exp(1500(T e - T g )/{T e T g )). 



L x = 4 mm and height L y = 4 mm. It consists of conductive copper strips (with negligible thickness) of width w = 1 
mm, separated by a 0.065 cm thick dielectric with width equal to 3 mm and relative dielectric permittivity e r = 5. 
The electrical capacity of the reactor is given by the conventional formula C — e r e S/d. 



B. Transport parameters and rate coefficients 

The working gas is a "airlike" mixture of a fixed fraction of nitrogen ([N 2 ]/N= .78) and oxygen ([0 2 ]/N= 0.22), 
as is normally present at sea level at p = latm. The electron homogeneous Boltzmann equation is solved with the 
2-term expansion in spherical harmonics [29] for a mixture of N 2 — 2 (22%), assumed thoroughly constant. The gas 
temperature is also assumed constant, both spatially and in time frame, with T g = 300 K, and the same applies to 
the vibrational temperature of nitrogen T V (N 2 ) — 2000 K and oxygen T v (0 2 ) = 2000 K. This assumption avoids the 
need of a more complex vibrational kinetic model. 

Using the set of cross sections of excitation by electron impact taken from Ref. [30j rates coefficients and transport 
parameters needed for the electronic kinetics arc obtained. So far, the species included in the present model are 
the following: N 2 , 2 , 2 and electrons. Remark that at 1 atm the concentrations of N± ions need to be 

introduced since they are bigger than , due partially to the reaction N% + N 2 — > N± , which occurs at a higher 
rate than the direct ionization. Also, notice that at p < 10~ 2 Torr loss by ambipolar diffusion becomes overwhelming 
important. Ion diffusion and mobility coefficients were taken from the Report [3l|: fj. Q -.N = 6.85 x 10 21 V -1 m _1 

s" 1 (on the range of E/N with interest here), ^ Q +N = 6.91 x 10 21 V -1 m -1 s^ 1 , and fi N +N = 5.37 x 10 21 V" 1 
m _1 s _1 . The gas density at p = 1 atm and assuming T g = 300 K is N — 2.447 x 10 25 #/m 3 . 

At atmospheric pressure the local equilibrium assumption holds and the transport coefficients {v^, v io 2 n , /x e , (i p , 
D e , Dp) depend on space and time (r,t) only through the local value of the electric field E(r, t); this is the so called 
hydrodynamic regime. 

To avoid the use of Navier-Stokes equations and to obtain a faster numerical solution of the present hydrodynamic 
problem it is assumed that the gas flow does not alter the plasma characteristics and is much smaller than the charged 
particle drift velocity. This assumption allows a simplified description of the flow. 

With the above assumptions charged species can be described by continuity equations and momentum transport 
equations in the drift-diffusion approximation. This last approximation is valid if their drift energy is negligible with 
respect to thermal energy. Also, in the drift-diffusion equation it is neglected the temperature gradient term (e.g., 
Ref. 0). 

The reactions included in the present "effective" kinetic model are listed in Table HI It is assumed that all volume 
ionization is due to electron-impact ionization from the ground state and the kinetic set consists basically in ionization, 
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FIG. 2: Program flow chart of the EHD code (CODEHD). 



attachment and recombination processes. The kinetics of excited states and heavy neutral species is not considered, 
particularly the possible important role of metastable species (certainly improving plasma density through Penning 
ionization (32|) and nitrous oxide (poisoning gas) are neglected. The major drawback associated to the present 
"effective" kinetic model is the lack of consistency due to the neglect of atomic species, but it allows to investigate 
the general trends of the plasma actuator operating characteristics. 

For more concise notation, we put n p2 = [N^ ]; n p \ = [N^ ]; n p = [0 2 ]; n n = [0%], and n e = [e]. 
With the above assumptions charged species can be described by continuity equations and momentum transport 
equations in the drift-diffusion approximation. This last approximation is valid if their drift energy is negligible with 
respect to thermal energy. Also, in the drift-diffusion equation it is neglected the temperature gradient term. 



C. Numerical model 



The particle's governing equations are of convection-diffusion type and they are solved using a well-known finite- 
element method proposed by Patankar [33|. Poisson equation for the electric field is solved applying the successive 
over-relaxation method (SOR) for boundary-value problem. The chosen time step is limited by the value of the 
dielectric relaxation time. For the present calculations it was used (100x100) computational meshes. The entire set of 
equations are integrated successively in time supposing the electric field constant during each time step, obtaining a 
new value of the electric field after the end of each time step. The method used to integrate the continuity equations 
and Poisson equation was assured to be numerically stable, constraining the time step width to the well known 
Courant-Levy-Friedrich stability criterion. 

The governing equation for N^" (at atmospheric pressure the nitrogen ion predominant is N4") is: 

+ V • (n p2 v p2 ) = S% 2 N 2 n p iK icl - K ic2 S N2 Nn p2 - K r2 n p2 n e . (1) 
The governing equation for (one of the most mobile charged specie in the plasma) is: 

—gj- + V • {n pl v p2 ) = n e v^ 2 n + K lc2 5 N2 Nn p2 - Pun n n pl - (3n e n pl - K icl S% 2 N 2 n pl . (2) 
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FIG. 3: Electric current, applied voltage, gas voltage and memory voltage as a function of time. Conditions: Case I. Solid 
curve: current; dot curve: V m ; dashed-dot curve: V 9 ; dashed curve: V . Conditions: V r ms=5 kV, f=5 kHz, dielectric width=3 



The oxygen positive ion considered is and its resultant governing equation is given by: 



Of 



f + V ■ (n p Vp) = n e v io 2 n - Pun n n p - f3n e n p . 



(3) 



The negative ion 2 was introduced (due to high electronegativity of oxygen), and its governing equation was written 
as: 



+ V • (n„v„) = i>au n e ~ Pun p in n - K da So 2 Nn n . 
Finally, the governing equation for electrons can be written in the form: 



(4) 



+ V ■ (n e v e ) = n e (v io 2 n + v io 2 n - vj t ) - /3n e (n p + n pl ) + K da So 2 Nn n - K r2 n p in e 



(5) 



Photo- ionization was not included as a non-local secondary effect. 

In order to avoid the use of Navier-Stokes equation, the drift- diffusion approximation for the charged particle mean 
velocities appearing in the continuity equations is used instead: 



riiVi = n^iE - V(niDi). 



(6) 



Here, fii and Di represent the charged particle mobility and the respective diffusion coefficient. The applied voltage 
has a sinusoidal wave form: 



V(t) = V dc + VbsimV), 



(7) 



where V dc is the dc bias voltage (although here we fixed to ground, V dc — 0) and id is the applied angular frequency. 
Vb is the maximum amplitude with the root mean square voltage in this case of study V rms = 5 kV and the applied 
frequency / = 5 kHz. 

In an RF discharge the calculation of the displacement current is a not a simple matter. However, the total current 
(convective plus displacement current) was determined using the following equation given by Sato and Murray [34j : 



FIG. 4: Electric field distribution in the X-Y plane at a given instant of time during the first-half cycle. 



t / x e f ( ^ 9n„ dn e dn n \ e f ( <9El _. \ , /0 , 

W) = vj v (np^P ~ neWe - n " w ™ - D p~q7 + e ~£h + n Jh ) ' L + VJ V (~dT ' L ) ' ( ' 

where J v dv is the volume occupied by the discharge, is the space-charge free component of the electric held. The 
last integral when applied to our geometry gives the displacement current component 

I dlsp {t) = - 2 — j v dv. (9) 

We assume throughout the calculations 7 = 5 x 10~ 2 since Auger electrons are assumed to be produced by impact 
of positive ions on the cathode with such an efficiency. Hence, the flux density of secondary electrons out of the 
cathode is given by 

3se(t) = 7j P (t), (10) 

with jj, denoting the flux density of positive ions. Secondary electron emission plays a fundamental role on the working 
of the asymmetric DBD plasma actuator. The progressive accumulation of electric charges over the dielectric surface 
develops a so-called "memory voltage" , whose expression is given by: 

1 r* 

V m (t) = JT- / I d (t')dt' + V m (t ). (11) 

Here, Cd s is the equivalent capacitance of the discharge. 

As charged particles are produced in the plasma volume, it is calculated the space-charge electric field by solving 
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the Poisson equation coupled the particle's governing equations: 

AV = - — (n p -n e -n n ). (12) 

We assumed the following boundary conditions: 

• over the electrode (Dirichlet boundary condition): V(x, y = 0, t) = V — V m \ 

• over the insulator (Neumann boundary condition): E n = (E • n) = 

The flux of electric charges impinging on the dielectric surface builds up a surface charge density a which was calculated 
by balancing the flux to the dielectric and it is governed by 

^ = e (|r p ,„|-|r e ,„|). (13) 

Here, r PjI1 and r e ra represent the normal component of the flux of positive and negative ions and electrons to the 
dielectric surface. Furthermore, it is assumed that ions and electrons recombine instantaneously on the perfectly 
absorbing surface. As it will be discussed later, this simplified assumption constitute a drawback of the present 
model. 

The entire set of equations were solved together self-consistently at each 1 ns time step, as illustrated with the 
program flow chart shown in Fig. [51 

III. RESULTS 

The electrohydrodynamic area of research has grown to a large extent lately, but it still remains to achieve a better 
understanding on how the charged particles transfer momentum to neutrals, and what physical limitations restrain 
the applicability of this device for boundary control and neutral flow propulsion. 

The simulations were done for a two-dimensional flat staggered geometry, while assuming the plasma homogeneous 
along the OZ-axis (see Fig. [T]). This is essentially a "surface discharge" arrangement with asymmetric electrodes. 

A. Electrical characteristics 

Fig. [3] shows the evolution along a full period of the calculated electric current (convective plus displacement 
currents), applied voltage, gas voltage and memory voltage. A stationary solution of the entire set of equations 
modeling the apparatus is achieved typically after 4-5 cycles with CODEHD. With the assumed conditions at about 
700 Volts (or V g — 2023 Volts) electron avalanches develop (the first avalanches occur in the gas gap), replenishing 
the volume above the surface with charged particles. In a coplanar discharge, the current density increases with the 
applied voltage and almost nearly with the permittivity [35| . 

This corresponds to roughly Ei on = V g /d = 31 kV/cm, when the theoretical dielectric breakdown of air is 30 
kV/cm. Hence, the charged particles flowing to the dielectric surface start accumulating on the surface, building-up 
an electric field that prevents the occurrence of a high current, and therefore quenching the discharge development 
at an early stage. At about 3 fis develops the first avalanche. The pulse form is typical: the first peak corresponds 
to the arrival of fast electrons to the electrode and it corresponds to the instant when the electric field is maximal; 
afterwards, the pulse slowly fade out, due to the dielectric charging and consequent decrease of the electric field. 

Also, from Fig. [5] it is clear that the electric field has a maximum peak at the edge of the exposed (and energized) 
electrode. From this edge the successive streamers (current pulses) are initiated. 

At 4th cycle the transient regime is totally damped and the electrical characteristics stabilize. During the first 
(positive) half-cycle electrons and negative ions are formed above the dielectric surface, which acts as a virtual 
cathode (buried below the dielectric), while they are pushed by the electric field towards the positive (energized) 
electrode located on the top of the dielectric playing the role of the anode, as can be seen from the set of snapshots 
shown from Figs. [5] and [71 The dielectric surface is progressively charged. 

Fig. [5] shows the surface charge density vs. time, for the case of a larger dielectric width. The charge density should 
be of the order a ~ 10~ 8 C/m 2 , since Q = CV — <rS, or a = CV/S = e e r V/d. It becomes increasingly negative 
during the first half-cycle and tends to reverse sign during the second half-cycle. The memory voltage is regulated by 
the surface charge density deposited over the dielectric by the successive streamers. Notice that when the energized 
electrode is negative, more streamers occur. 




Time (us) 



FIG. 5: Surface charge density as a function of time. Conditions: Vrms=5 kV, f=5 kHz, dielectric width=3 mm. 




FIG. 6: A time series shown of the electrons density at the threshold of an avalanche during the first half-cycle. 



FIG. 7: Two dimensional distribution of N% ions particle density predicted by the numerical code for typical conditions, during 
an avalanche occurring at the first half-cycle. 

Fig. [6] shows a snapshot of the electrons density along time since the beginning of an electron streamer occurring 
during the first-half cycle. The density profiles were obtained after 4-5 cycles with the dielectric charging onset 
resulting in the formation of a virtual negative electrode. Electrons are driven along the electric field lines toward the 
anode (the energized electrode at the first-half cycle). 

At the same time and in the reversed direction a streamer of positive ions N£ are driven to the cathode, as it 
is shown at the time series illustrated in Fig. [7J During the streamers ponderomotive forces attain their highest 
magnitude. The event occurring at the instant of time t — 6.75/xs portrays the arrival of ions to the virtual 
cathode. 

We notice that the ions are mainly created in volume and are immediately driven to the virtual cathode, while 
an average number of them are retained nearby the anode decreasing almost linearly with distance from the energized 
electrode surface. Thus, it is clear that the propagation of the electron avalanche near (and above) the dielectric 
surface is of considerable importance, dictating the strength of the body forces a nd g as speed. 



Conditions for maximizing ion-driven gas flow were obtained by Rickard et at [361 ]. and they concluded that, irre- 
spective of geometry, ponderomotive forces on the gas are maximized by increasing current density and by decreasing 
mobility (i.e., charge carriers which exercise highest drag on the neutral gas). Therefore, A^ 1 " ions seem to be the best 
candidate for this purpose. 

Our numerical model shows that when the present conditions prevail, heavy species such as N^" move more slowly 
with the varying external electric field. In fact, the most actives species on the process of momentum transfer are the 
electrons and N^~, but the molecular ions due to their mass contribute in majority to control the boundary layer and 
propulsive force. 

From Fig. [JJ it is noticeable the multiplications of ions when flowing from the electrode edge to the dielectric 
surface, flowing along the reverse way as electrons did. The ions feeding along the dielectric surface is due to a relative 
bigger dielectric width which favors the increase of the ions swarm [3 71 ] , and thus increasing the gas speed due to the 
momentum exchange onset from charged particles to neutrals. Notice that at t = 6/Lts nitrogen ions leave the region 
at the boundary between the electrode and the dielectric, which corresponds a region of maximum electric field (see 
Fig. E} . 

Fig.[8]shows contours of constant potential at t = 5/xs. We can see the decrease of the potential above the energized 
electrode and a field reversal region toward the dielectric side (with the negative electrode below). The negative glow 
remains in the proximity of this region, while a second region of field reversal is also momentously observed. The 
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FIG. 8: Contours of constant potential at t = 5/xs. Same conditions as in Fig. 3. Compare with Figs. 4-5. 



region of negative electric potential that appears at 5 /xs in the first half-cycle is due to the presence of an excess of 
negative charge due to ions, while the region of maximal electric potential has as field source all others positive 
ions. Numerical simulations have shown that while negative ions in the air do not contribute significantly to the 
ponderomotive forces, they can play a role in the discharge working (25j . Another major structure of the normal 
glow discharge remains above with higher values of the potential. Hence, the phenomenology and typical structures 
developed by normal glow discharge are also displayed by the OAUGDP™. These aspects were also shown in previous 
publications, like the one-dimensional numerical simulations of the OAUGDP™ done by Ben Gadri [2?], HH and fast 
photography obtained by Massines et al. [39] . 

Hence, the field reversals observed have the strong contribution of the negative ions which are thus playing an 
important role in the discharge working, since it is well known that field reversal location decides the fraction of ions 
to the cathode and the magnitude of the plasma maximum density [H, l40l [ill |42| . 

The EHD force acting on the charged particles is given by [18|, [22[ : 

F = e(n t - n e )E - V(mkTi + n e kT e ) + (m e u e - rrnu^S, (14) 

where rii and n e are respectively the ion and electron number density; u e and are the electron and ion mean 
velocities; S is the charged particle production rate. The last term of Eq. [TJ] is here neglected, since its importance 
confines to other phenomena such as electrophoresis and cathophoresis 43)1 . and its contribution to the induced 
velocity is smaller due to the very small electron/neutral atoms mass ratio 44|. Although we had assumed a constant 
ion temperature in equilibrium with the gas temperature, we calculated the electron temperature and verified that 
the order of magnitude of the second term is about 1 % of the coulombian force term which constitute the main 
term of the theory of paraelectric gas flow control developed by Roth [19(. We calculated the electrostriction force 
term (not included in Eq. [T4"]) and concluded that is not significant, contributing at maximum with 1 % to the total 
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FIG. 9: Electrohydrodynamic forces (in N/m) as a function of time. Same conditions as in Fig. [3] 



ponderomotive force. Subsequently, the ponderomotive forces were averaged over the area of calculation. 

It is found that the calculated space averaged ponderomotive forces per unit volume increases when the electrode 
width increases @. On average, during the second half-cycle the ponderomotive force magnitude decreases with a 
magnitude of a few N/m as shown in Fig. [21 a result consistent with experimental results as such presented in Ref. [45| . 
This happens when the voltage polarity is reversed and the energized electrode play the role of cathode. This is due 
to the potential gradient reduction on the edge of the expanding plasma (see also Ref. Q). This numerical result is 
contrary to experimental study presented in Ref. [461 ] . showing that the forces decay exponentially with increasing 
electrode diameter. This is due to the role of the dielectric that in our model was assumed to absorb electric charged 
particles instead to feed the swarms and strengthening the body forces (see also Ref. B). However, there is still 
no consensus on the ponderomotive force dependency. For example, Singh and Roy [47| obtained the magnitude 
of approximated force and have shown that it increases with the fourth power of the amplitude of the rf potential, 
implying that the induced fluid velocity also increase. This is certainly an aspect that must be dealt with more 
caution. 

In fact, electrons are faster than ions. After the first breakdown, they start to charge negatively the dielectric 
during the 1st half-cycle; positive ions gains more energy and electrons are also increasingly accelerated, due to a 
higher growing potential. Therefore, the positive ions density is bigger during this half-cycle, resulting in stronger 
ponderomotive forces. 

During the 2nd half-cycle, positive ions (which are mainly formed in volume) tend towards both the electrode and 
the dielectric surface; the charge surface density on the dielectric start to become less negative and the gas voltage 
decrease, generating less ions and electrons, resulting in smaller body forces; this is the mechanism of an asymmetric 
flat panel device. 

Otherwise, if the discharge is entirely symmetrical in both half-cycles, it is expected that the average gas speed 
equals zero. In fact, it is the asymmetry in the streamers (in our case bigger during the fist half-cycle) that gives an 
overall positive gas speed along the axis. 

In our present model, calculations of EHD ponderomotive force have shown that its maximum intensity is attained 
during electron avalanches, with typical values on the order of 5 x 10 9 N/m 3 . F x points along OX (propelling direction), 
while F y points downwards (boundary layer control). Our calculations show that the resulting average gas speed is 
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FIG. 10: Gas speed (in m/s) as a function of time. Same conditions as in Fig. [9] 

about 20 m/s and the net EHD body forces (with the present conditions) have comparable values in the first and 
second-half cycle, although slightly bigger during the first half-cycle, as shown in Fig. |T0l It is clear that the successive 
streamers that charge the dielectric surface are responsible for pulling the flow upstream, unidirectionally, as recent 
experiments have shown [45[. 

IV. CONCLUSION 

A two-dimensional fluid model of an asymmetric plasma actuator display the behavior of charged species during 
both half-cycles when electrodes are subject to a sinusoidal applied voltage. The actuator is strongly dominated by 
dynamics, charged species form preferentially at the edge of the electrode with the insulator, and their subsequent 
behavior and ability to provide an unidirectional gas speed results from the interaction of the charged species with the 
dielectric, in particular, the effect of the electric field above the insulator and the propensity of the dielectric surface 
to adsorb or not charged species, and thus controlling the plasma density of the streamers. An appropriate model to 
describe a realistic interactions of charged species with the dielectric for plasma density enhancement remains to be 
done, lacking in the literature a more careful study of this important issue. 
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